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1. Introduction 

In its minimal version, the Standard Model of particle interactions requires a complex scalar 
weak-isospin doublet to spontaneously break the electro-weak gauge symmetry. Three of 
its four original degrees of freedom transform into the longitudinal modes of the and 
Z gauge bosons while the fourth manifests itself as a neutral scalar field, the Higgs boson. 
Yukawa-type couplings between the Higgs field and the fermions generate gauge invariant 
masses for the fermions. 

There is no strong theoretical reason to believe that the minimal version of the Stan- 
dard Model is indeed realized in nature. Extended models give rise to several Higgs parti- 
cles which can be electrically charged or neutral and which can be even or odd under CP 
inversion. 

The most appealing extension of the Standard Model, at the moment, is the Minimal 
Supersymmetric Standard Model (MSSM), which basically doubles the field content of the 
Standard Model and requires two Higgs doublets for the generation of the particle masses. 
This results in five physical Higgs bosons: three neutral {h,H,A) and two charged (H^). 
h and H are CP-even while A is CP-odd (for a review see Ref. |j^]). 

To date, no Higgs boson has been observed in nature, in spite of great efforts that 
have been made at particle accelerators, especially at the Large Electron Positron collider 
(LEP) at CERN. The null results lead to lower limits on the mass of possible Higgs bosons. 
For a scalar Higgs boson in the framework of the minimal Standard Model, this limit is 
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Mh > 114.4 GeV at 95% CL. Assuming the MSSM, the hmit for a neutral scalar Higgs goes 
down to about 91 GeV. The 95% CL limit on a CP-odd Higgs boson, the subject of this 
paper, is around 92 GeV ||2|. 

The Large Hadron Collider (LHC) at CERN is scheduled to commence taking data in the 
year 2007. It will be a proton-proton collider with a center of mass energy of ^/s = 14 TeV 
and has been designed specifically to search for Higgs bosons. For all relevant values of 
the Higgs boson mass and most kinds of Higgs bosons in various models, the gluon-gluon 
production mode will be of the greatest importance both for discovery and for measuring 
the Higgs boson mass. For CP-even Higgs bosons, the theoretical prediction is now fairly 
well under control, having been calculated to next-to-next-to-leading order (NNLO) in the 
strong coupling constant ^. 

In this paper, we use the techniques of Ref. [^, ^, |6| to evaluate the production rate for 
a pseudo-scalar Higgs boson to the same accuracy, i.e., to next-to-next-to- leading order in 
QCD. We work in the heavy-top limit, using an effective lagrangian for the interaction of 
the pseudo-scalar Higgs boson with the gluons. As in the case of the scalar Higgs boson, 
this does not restrict the validity of the results to the Higgs-mass region well below 2Mt if 
one factors in the full top mass dependence at leading order 

In the effective lagrangian approach, the massive top-quark loop that mediates the 
coupling between the Higgs and the constituents of the initial state hadrons, reduces to 
effective vertices with known coefficient functions. What remains to be computed at NNLO 
are 2 — > 1 processes up to two loops, 2 —> 2 up to one loop, and 2 ^ 3 at tree level, where 
all internal particles are massless, and all external particles are taken on-shell {pf = for 
quarks and gluons, = M\ for the Higgs). 

We note that we do not consider contributions from virtual particles in extended 
theories in this paper. We also do not consider the effect of b quark loops at this order, 
which can be important for example in the MSSM, when the coupling is enhanced by a large 
value of tan [3. For light quarks, like the 6, the effective lagrangian cannot be formulated 
and one must perform a true three-loop calculation at this order. That calculation is 
beyond the current state of the art. 

NLO corrections to this process were evaluated in Ref. |0, ^ and found to be very similar 
in size to the NLO effects for scalar Higgs production. We find that this is also true at NNLO: 
The K-factors at NNLO for the scalar and the pseudo-scalar are comparable. This means 
in turn that the production rate for pseudo-scalar Higgs bosons is fairly well under control. 
Though still sizable, the NNLO term in the perturbative series is significantly smaller than 
the NLO term and still higher order effects are presumably negligible. Accordingly, the 
unphysical dependence of the cross section on the renormalization and factorization scales 
is reasonably small, allowing for a prediction of the total rate with errors of the order of, 
or below, the expected experimental precision. 

The paper is organized as follows: In the second section, we describe our theoretical 
framework, including the effective lagrangian, its Wilson coefficients and the renormaliza- 
tion of its operators, and our prescription for handling Levi-Civita tensors and 75. In the 
third section, we briefly discuss the matrix element calculations and the calculation of the 
partonic cross sections. In section 4, we present our results for the partonic cross sections. 
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The result for pseudo-scalar production is very similar to that for scalar production and 
the expression for the difference between the two is quite compact. The full expression is 
presented in the appendix. In section 5 we compute the hadronic cross section by folding 
in the parton distributions and finally we present our conclusions. 



2. Theoretical setup 

In principle, the whole theoretical background is laid out very clearly in Ref. ||9|. Never- 
theless, let us review the necessary ingredients for our calculation. 

We assume that the pseudo-scalar Higgs boson. A, couples only to top quarks, t. The 
interaction vertex is given by 

A 

Cau = -igt- Mti-f5t , (2.1) 

V 

where gt is a coupling constant that depends on the specific theory under consideration. 
In the MSSM, for example, one has gt = cot /3, where tan /3 is the standard ratio of the 
vacuum expectation values of the two Higgs fields. 

The coupling of A to two gluons is then mediated by a top quark loop. In the heavy-top 
limit, this interaction can be described by an effective lagrangian M: 



^Agg = -gt — C'f of + C2 ^2 



V 



of = G%G^'^^ , Of = 5 J 5] 759 



(2.2) 



G^j^ is the gluon field strength tensor, and G^^, is its dual: 

^aIi/ — ^fJ.val3 G""'"^ . (2-3) 

All quantities in these equations are to be understood as bare quantities in the effective 
theory of five massless flavors. Thus the sum over quarks (X]q) ™ Eq. ( |2.2| ) does not include 
the top quark. Cf and Cf are coefficient functions that can be evaluated perturbatively. 
One may define renormalized operators and coefficient functions as follows: 

2 2 
6, = ^Z,pf, a = Y,iZ-%Cf, i = l,2. (2.4) 



The renormalization matrix Z is given by pA\ 



Z21 = , Z22 = ^MS -^1 ) 
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with 
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where ?^/ is the number of hght (i.e. massless in our approach) quark flavors. In our numer- 
ical analysis below, we will always assume nj = 5. Zq,^ and Z^g are the renormalization 
constants of the strong coupling and the singlet axial current, respectively. Z| is a finite 
renormalization constant that will be discussed below. 

The coefficient functions are process independent and have been evaluated in the con- 
text of pseudo-scalar Higgs decay to NNLO M: 



' vr 16 



a.o 



vr 



1 , 

— in — K 

4 M2 



+ 0(a3) 



(2.7) 



where Mt is the pole mass of the top quark. 

The matrix elements to be evaluated have the form 



(CiOi + 0262) = CiZu {Of) + {CiZu + C2Z22) {Of) , (2.8) 

with {O^) = {ab\0^\XH), where a and b label the two partons in the initial state, and X 
denotes an arbitrary number of partons in the final state. We need to evaluate the square 
of this matrix element up to 0(af). 






(a) 




id) 



Figure 1: Sample diagrams for (Of). Higher orders are obtained by dressing these diagrams with 
additional quarks and gluons, both virtual and real. 



of generates vertices which couple two or three gluons to the pseudo-scalar Higgs 
(the ggggA-veitex vanishes due to the Jacobi-identity of the structure functions of SU(3)). 
Thus, the diagrams related to {Of) are the same as in the scalar Higgs case (cf . Ref . [|^] ) . 
A sample of typical diagrams is shown in Fig.|l|. Higher order contributions are obtained 
by dressing these diagrams with additional gluons and quarks, either as virtual particles, or 
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(a) (b) (c) 

Figure 2: Sample diagrams for (Of) as required up to NNLO 



^^^^ 

id) 



as real particles in the final state. The actual results for these diagrams are different from 
the ones in the scalar case, of course, due to the different Feynman rules corresponding to 
G^,G^''' and G^.G'^^ 

As for (of), since its prefactor in Eq. (|2.8D is of order a^, it is only required to one 
order less than (Of). Diagrams contributing to (Of) for the different subprocesses are 
shown in Fig. ^. They appear in interference terms with the corresponding diagrams of 
{Of) (see Fig. I). At NNLO, the square of diagram (c) is the only potential contribution 
to the process arising solely through O^. It vanishes, however, in the limit that the light 
quarks are massless. 

When evaluating the diagrams, the presence of the manifestly 4-dimensional Levi- 
Civita tensor in Of as well as the 75 in O^ requires special care. We strictly follow the 



strategy outlined in Ref. [11|. This means that first we replace 



3! 



(2S 



Since we need to evaluate squared amplitudes, there will always be exactly two Levi-Civita 
tensors in our expressions. This allows us to express them in terms of metric tensors, which 
can be interpreted as d-dimensional objects: 



(2.10) 



where, as indicated, the brackets around the upper indices indicate that their positions are 
to be anti-symmetrized while the lower indices remain fixed. All further calculations (i.e., 
integrations, contraction of indices and spinor algebra) can be performed in d dimensions. 
This procedure is slightly different from that of Ref. , but we have checked that it gives 
the same result for the pseudo-scalar Higgs decay up to 0(a^).^ Furthermore, we have 
computed the two-loop virtual corrections for our process in both approaches and obtained 
identical results. 

The finite renormalization constant Z|, introduced in Eq. ( |2.6D , is determined by re- 
quiring that the one-loop character of the operator relation of the axial anomaly is preserved 
also at higher orders: 



TT O 

Note that 0\ and O2 denote renormalized operators here, as defined in Eq. (I 



(2.11) 



^We thank M. Steinhauser for providing us with unpublished intermediate results concerning Ref. 
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3. Methods of evaluation 



3.1 Virtual corrections 

Potential sub-processes are gg ^ A and qq — > A. It turns out, however, that the latter 
does not contribute at NNLO. The diagrams needed are two- loop vertices with one massive 
and two massless external legs. 



For their evaluation, we use the technique of Ref. |12] that has already been applied 
successfully to the evaluation of the virtual correction to scalar Higgs production Q . This 
means that the two-loop vertex diagrams are first mapped onto three-loop propagator dia- 
grams by interpreting the massless external lines (gluons or quarks) as part of an additional 
loop. The resulting integrals, including tensor structures, can be treated by means of the 



integration-by-parts algorithm of Ref. [y^]. In particular, we can use the FORM [14| program 
MINCER as the basis of an algebraic program that reduces all integrals encountered 
to a set of master integrals [Sj. The analytic expressions for the latter have been known 



for a long time |16|. Note that a generalized version of the method of Ref. [|l^ ] has been 



constructed in Ref. 

3.2 Single real emission 

The radiation of one additional parton has to be evaluated up to one-loop level. The 
contributing processes are gg — > Ag, gq Aq, gq — > Aq, and qq Ag. The one- 
loop matrix elements have been evaluated to all orders in the dimensional regularization 
parameter e = (4— (i)/2. After interfering with the tree-level amplitude, the squared matrix 
element can be integrated over single-emission phase space to obtain the contribution to 
the partonic cross section in closed analytic form. The interference of bare operators Of 
and of is of order e. Nonetheless, operator Of contributes to the single real emission 
cross section through operator mixing since Z12 is of order 1/e. The Feynman rules, loop 
integrals and phase space integration have all been implemented in FORM programs. 



3.3 Double real emission 

We need the tree-level expressions for the processes gg — > Agg, gg Aqq, gq — > Agq, 
qq — > Agg, qq — > Aqq, and qq Aqq (and the corresponding charged conjugated processes). 
The squared matrix elements can be evaluated straightforwardly in d = 4 — 2e space-time 
dimensions with the help of FORM. The result is a rather large expression of several thousand 
terms which must then be integrated over phase space. 

The phase space for double real emission is quite complicated and we perform the 
integration in the method of Ref. Q. That is, the matrix element and the phase space are 
expanded in terms of (1 — x), where x = M\/s. The result is a power series expansion in 
(1 — x) and ln(l — x) (at NNLO, the highest power of the logarithm is ln'^(l — x)) for the 
double real emission contribution to the hadronic cross section. If one were to compute all 
terms in the series, this would be an exact result. 

In fact, a truncated series is sufficient for obtaining the NNLO cross section to very 
high numerical precision. In this case, one obtains a cancellation of infrared singularities 
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by expanding the other contributions to the partonic cross section (single real emission, 
mass factorization; the dependence of the virtual terms on x is simply oc (5(1 — x)) to the 
same order in (1 — x) as the double real emission term. In Ref. we computed scalar 
Higgs boson production to order (1 — x)^^ and arrived at a prediction for the hadronic cross 
section that is phenomenologically equivalent to one based on the closed analytic form of 
the partonic cross section that has recently become available jQ]. This comes as no great 
surprise. The functions which contribute to the closed analytic form can all be expanded 
in terms of (1 — x). For example, the dilogarithm can be represented as 

^-^ n \n I 

n=l ^ ' 

Furthermore, steeply falling parton distributions ensure that the threshold region domi- 
nates and that convergence in (1 — x) is quite rapid. 

Although the truncated series leads to a physical result that is by all means equivalent 
to the exact expression, the approach of expanding in (1 — x) can be taken farther. If the 
expansion can be evaluated up to sufficiently high order in (1 — x), one can actually invert 
the series and obtain the partonic cross section in closed analytic form. This is due to the 
fact that only a limited number of functions appear in this closed form representation: log- 
arithms, dilogarithms, and trilogarithms of various arguments, multiplied by 1/x, 1/(1 + x) 
and (1 — x)",!! = —1,0,1,2,3. Taking the functions which appear in the result for the 
Drell-Yan cross section and these possible factors, one finds that carrying out the ex- 
pansion to order (1 — x)^^ should suffice to allow inversion of the series. In fact, we have 
carried out the expansion to order (1 — x)^*^*^ so that we could over-determine the system. 
As a check, this procedure has also been carried out for the scalar Higgs boson Q, and 
complete agreement with the result of Ref. Q was found. 

It is interesting to observe that the task of inverting the series is much simpler when 
one examines only the difference between the scalar and pseudo-scalar cases since they are 
already very similar at the level of squared amplitudes. The difference may be inverted 
with less than twenty terms. 

4. Partonic results 

As noted above, the partonic cross sections for scalar and pseudo-scalar Higgs production 
are very similar, so that many terms cancel in the difference between the two. Thus, the 
partonic cross section for pseudo-scalar Higgs production can be expressed conveniently in 
terms of the known scalar Higgs boson cross section (with modified normalization) plus a 
remainder. For this purpose, we write 

A„fe$, «>G{F,^}, a,h^{g,q,q], (4.1) 

where dah^ is the cross section for the process a6 — > $ + X. o and h label the partons in 
the initial state, $ means either a scalar (i?) or pseudo-scalar (A) Higgs boson, and X 
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denotes any number of quarks or gluons in the final state. In the normalization factors, we 
keep the full top mass dependence: 



TT 



TH [1 + (1 - Th) f{TH)\ 
2 

gt TA f{TA) , r$ 



4M2 

a4 ' 



(4.2) 



where Gp ~ 1.664 x 10 ^ GeV ^, and gt has been introduced in Eq. (2.1). The one-loop 
function appearing in Eq. (J^^) is defined by 



fir) 



arcsm 



2 1 



In 



ITT 



r>l, 
, r < 1. 



(4.3) 



For completeness, we note that in the heavy-top limit, the normalization factors approach 
the values 



Mt^p o TT a/2 Gf 

' 576 



0"A 



Mt^oo 



ttViGf 
256 



\9t\ 



(4.4) 



The kinematic terms are written as a perturbative expansion: 



Kb^ix) = 5ag5,g 5{l - x) + ^ A« (x) + (^) a2!,(x) + 0{ 



TT 



\ TT 



a. 



(4.5) 



The results for scalar Higgs production up to NNLO = H) can be found in Refs. |p^ , 
ii- 

With the help of these expressions, the corresponding results for the pseudo-scalar 
Higgs production cross section can be written in a rather compact form. At NLO, the 
result has been evaluated some time ago : 



^11!4(^) = ^IIh(^) + \ ^aghg 5{l - x) , a,b e {g, q, q} 



(4.6) 
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The main result of this paper are the NNLO terms in IS.ahA- We find: 



(2) 



1939 19 



144 



5(1 -x) + 6 Pi (x) 



(12 X - 6x2 + 6x^) ^^^^ _ _ 9a; ln2(x) 



3 (10 - X - 13x^ + 4x3 
^ 2 1 -X 



2x^) , , , (154 - 189x + 24x2 + 11x3) 
ln(x) H 



+ nf 



13 2 \ 2 

T- - - + 2^2 (5(1 - x) + - X ln2(x) + X ln(x) 

lb 3 / 3 



(1 - llx + 10x2 



6 



with 



1 1 , 
"4 + 2^^*' 



A 



(2) 
gqH 



, , (4-4x + 2x2) , 28 , 

^x) H in(l — x) — — X m (x) 

3 9 



(4.7) 



+ 



(22 + 30x-x2) (337 - 382 x + 51x2) 



■ ln(x) + 



18 



.(2) 



— X ln2(x) + — (6 + 11 x) ln(x) + — (37 - 40 x + 3 x2) 



(2) , 32 ^ ,2/ ^ . 32 



A,^a(^) = A,?h(^) + ^ ^ (x) + - (3 + 8 x) ln(x) 



a(5^(x) = A;;^,^(x) - ^ X ln2(x) + ^ (2 + 3x) ln(x) + ^ (11 - 12x + x2) , 

where Im = ln(M|/Mj?), ^„ = ({n) is Riemann's C function (C2 = n^/G = 1.64483..., 
Cs = 1.20206...) and2?n(x) = [ln"(l-x)/(l-x)]+. Of course one must change Mjy ^ Ma 
wherever it appears in AabH- For the sake of completeness, we list the full result for AabA 
in App.^. The difference A^bA ~ AabH can be expanded readily in terms of (1 — x) in 
order to bring it to a form consistent with the results of Ref . Q . 

The only contribution that originates from the presence of the (renormalized!) operator 
O2 is the term 82 in the equations above. We obtain it by computing the interference of 
diagram (a) in Fig. || with diagram (a) of Fig. ||. Alternatively, it can be derived from 
Eq. ( p. 11 ) by using the LO result for AggA- 

as uf /as^ 2 

X 



16 , 
+ ^(11 — 

.(2) 



27 

9x2 



X'^) + Uf 



— X m x) (1 

27 ^ ^ 27 ^ 



16 



16 



(4.8) 



vr 8 L " ^ ^ 

which leads to the same result for 62 as above. This provides a welcome check on the 
normalization of the contribution from ©2. 



5. Hadronic results 



In complete analogy to scalar Higgs production, one has to convolute the partonic rate 
with the proper parton distribution functions, in order to arrive at a physical prediction 
for the hadronic production rate: 

(^hih2{s)=Y^ dr -[fa/h^{Xa)fb/h2i'r/Xa)]^ab{s = ST). (5.I) 

a,b 
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A consistent NNLO hadronic result requires not only a NNLO partonic cross section, but 
also PDFs that have been evaluated at NNLO. Strictly speaking, such a set of PDFs is not 
yet available, because the NNLO evolution equation is not fully known. Nevertheless, we 
use an approximate set of PDFs based on approximations to the evolution equation |19| 
derived from the available moments of the structure functions [201. At lower order, we use 



the corresponding lower order PDFs [21 
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Figure 3: Total cross section for inclusive production of pseudo-scalar Higgs bosons (A) at (a) the 
LHC {^/s = 14TeV) and (b) Tevatron Run II = 2TcV). The coupling constant for the coupling 
of A to top quarks is such that gt = 1. The cross section for other values of gt (e.g., gt = cot/3 in 
the MSSM) can be obtained by scaling the curves with \gt\^. 
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Figure 4: Total cross section for inclusive production of pseudo-scalar Higgs bosons (A) at the 
LHC (v^ — 14TeV), as in Fig.^ for a larger mass range. The cusp is an effect of the top-quark 
threshold. 
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Figure 5: Scale dependence of the cross section for an Ma = 115 GeV pseudo-scalar Higgs boson 
at the LHC (y^ = 14TeV) (a) varying fj,p = /i^, (6) varying ^p, fJ-B. = Ma and (c) varying /i;?, 
= Ma. 
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Figure 6: Same as Fig. S for the Tevatron Run II (-/s = 2TeV). 



Fig. ^ shows the total cross section as a function of the pseudo-scalar Higgs boson mass, 
Ma, at (a) the LHC, and (6) Tevatron Run II. Fig. ^ shows the curves for the LHC over a 
larger mass range. The renormalization and factorization scales have been set to M4. The 
cusp in Fig. ^ is a leading order effect caused by the tt threshold, cf. Eq. (|4.2| ). Since we do 
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not want to specify a particular extension of the Standard Model, we choose the coupling 
of the pseudo-scalar Higgs to the top quarks such that gt = I (cf. Eq. ( |2.1| )). The total cross 
section scales with \gt\'^, so that the actual numbers can be easily obtained from the ones 
shown in the figure. As already noted, the behavior of the corrections is very similar to 
the ones for scalar Higgs production: The NNLO corrections are significantly smaller than 
the NLO corrections, indicating a nicely converging result with uncertainties well under 
control. 

This is further affirmed in Figs. ^ and |6|, which show the variation of the cross section 
with the renormalization and factorization scales for a fixed Higgs mass = 115 GeV. 
The scales are varied between M^/A and AM a'- in sub-panels (a), the renormalization 
scale and the factorization scale are identified and varied simultaneously; in sub-panels 
{h), the renormalization scale is identified with Ma^ and /ii? is varied, while in sub- 
panels (c), hf = Ma is fixed and fiR is varied. One observes a clear reduction of the scale 
dependence at NNLO with respect to NLO. Compared to the LO curve, there is a significant 
relative improvement in the scale variation, while the absolute ranges of variation at LO 
and NNLO are comparable. It is also clear that the scale dependence is dominated by 
renormalization scale dependence; factorization scale dependence is quite small. Using 
the rather conservative range of Ma/ A < fj, < AM a, one arrives at an uncertainty in a 
from scale variation of about ±30% (±65%) at LO, ±30% (±50%) NLO, and less than 
±20% (±25%) at NNLO for the LHC (numbers in brackets for Tevatron Run II). Varying ^ 
between Ma/2 and 2Ma results in a variation of a of ±20% (±40%) at LO, ±15% (±25%) 
at NLO, and less than ±10% (±15%) at NNLO. 



6. Conclusions 



The hadronic cross section for the production of a pseudo-scalar Higgs boson has been 
calculated at NNLO in QCD. We find that the corrections are similar to the scalar case, both 
in their magnitude and in their uncertainty due to scale dependence. While this uncertainty 
is still rather large, it seems that the NNLO calculation yields a reliable prediction for the 
total rate. 

Note Added: As we completed this manuscript, we became aware of a similar paper 
by Anastasiou and Melnikov|22]. We have compared results for the partonic cross sections 



and find complete agreement. 
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A. Analytic results 

In this appendix, the analytic expressions for the partonic cross section of pseudo-scalar 
Higgs production at NNLO are listed. Note that using the formulas of Eqs. ( [4.6[ ) and 
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( |4.7| ), one can transform these expressions into the corresponding ones for scalar Higgs 
production. 

At NLO 0, ^, the result for the gluon-gluon sub-process is 



A^g^^ix) = (6 + 6C2)<5(l - a;) + l2Vi{x) - (24x - 12 x'^ + 12x3) ln(l - x) 

, {l-x + xr , . , 11,, ,3 ^^-^^ 

— 6 ln(x) (1 — X) . 

1 — X 2 



For the quark-gluon channel one finds 



(1) (8-8x + 4x2) (4-4x + 2x2) (3-6x + x^) . 
^ota(^) = o - x) ln(x) , lA.2j 



and for the quark-anti-quark channel 



Ai;l4(x)=|(l-x)^ (A.3) 



The contributions from different sub-processes at NNLO are written as 



A^Z = ^a,A + nfA^,^, (A.4) 



where nj is the number of light (i.e., massless in our approach) quark flavors. 
For the sub-process with two gluons in the initial state, we find 
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741 139 _ 165 _ 9 ^ 

-g- + - - 2qC2 



5(1 - x) 



Po (x) + [139 - 90 C2] (x) - 33 V2 (x) + 72 P3 (x) 

- (144 X - 72x2 + 72 x^) ln=^(l - x) - (297 - 381 x + 348x2 - 330 x^) hi^{l - x) 

9 (31- 30x + 93x2 -94x3 + Six''), 

-^^ - x) ln(x) 

2 1 — X 

(2027 - 2735 X + 2182x2 - 2583x3) 



+ 



+ (180x-90x^ + 90x3)C2 



ln(l - x) 



„ (88 - 211 x + 312x2 -365x3 + 187x4) , 
+ 3 : m(l — x) m(x) 



+ 9 



+ 36 



- 18 



1-x 

(7 + 3x + 19x2 -3x3 -19x4 + 9x5) 
(1 - 6x -13x2- 6x3 + x4) 



In(l-x) ln2(x) 



1 + x 

(l + 2x + 3x2 + 2x3 + x4) 

1 + x 



ln(l-x)Li2(l-x) 
ln(l -x)Li2(l -x2) 



_ 9 (24 - 38x + 8x2 _^ 543,3 _ ;^9^4 .^9^5 
~ 2 1 -x2 

_ 9 (27 + 35x + 75x2 - 29x3 - 78x4 + 6x^ 
~ 2 l-x2 

9(l + 2x + 3x2+24x3 + 16x4)^. 2n 

9(l + 2x + 3x2-8x3-8x4) / (1-x 



-Li3(l - x) 



Lis 



(1-x) 



8 1 + x 

9 (7 + 14x + 21x2 + 8x3 + 4x4 
2 



Lis 



Lis 



x^ 

1 - X 



+ 
+ 
+ 



1 + x 

3 (317 - 398 X - 87 x2 + 300 x3 - 121 x4 

4 1-x 
9 (11 + 31x + 59x2 _ 25x3 - 65x4 ^ -^-^^5 
2 1 -x2 

(42 + 36x- 63x2 -33x3; 



1 + x 

^Li2(l-X 



- Lis - 



1-x 
1 + x 



-Li2(l — x) ln(x) 

~2 o 



^ , 2n 9(5 + 10x + 15x2-2x4)^. ^ 2^ 

Ll2(l-x2)+ ^T,_n 



4 ^ ' 4 1+x 

3 (21 + 23x + 41x2 - 37x3 -44x4 + 4 x5) ,^3 

4 1 - x^ 

3 (154 - 365 x + 675x2 - 827x3 + 374x4) , 

^ln^(x) 



-Li2(l-x2)ln(x) 



■ln3(x) 



1-x 



1 (2213 - 5599 x + 6603 x2 - 7003 x3 + 4342 x4) 



+ 9 



1 — X 

(9 - 2 X + 27x2 - 34 x3 + 9 x4) 



1 — X 



+ ^ (319 - 277x + 233x2 - 363x3)^2 - ^ (2x - x2 + x3)C3 



C2 ln(x) 

351 
2" 



ln(x) 



(16309 - 20611 X + 23819 x2 - 22749 x3) 



48 



(A.5) 
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and 



+ ■ 



689 , 5 . 5 . ' 
■^ + /i/t-3C2 + gC3 



(5(1 -x) 



2C2 



:Po(a;) 



10 



+ -(8 - 12 a; + 3 - 17x^) ln2(l - x) + -{x + x^) ln^{l - x) ln{x) 



(922 - 294x + 249 ar^ - 1570 
108 

2 (17 + 7x + 21 x2 - 61 a;3 + 25 



1 — X 



ln(l — x) 

ln(l — x) ln(x) 



16 



(x + x^) In^(x) ln(l — x) H (x + x^) ln(l — x) Li2(l — x) 

3 3 



+ 
+ 
+ 
+ 
+ 



(2 + 14x + 17x2)_ ^ (2-34x-31x2)^. 
Li3(l — x) L13 

1 (68 - 302x + 21x2 + 227x3 + 4x'i) 



(1 - x) 



X 



36 1-x 
(2-50x-47x2^ 



Li2(l-x) 



12 



Li2(l — x) ln(x) 



(2 + 6x + 9x2) 1 (68 + 100 x + 69x2 -351x3 + 132x4) ^ 



In'^fx) + 

72 ^ ^ 72 1-x 

1 (1282 - 382 X + 117x2 -3041x3 + 2384 x'' 



In^(x) 



216 1-x 
(12707 - 606x + 1641x2 - 17774x3) _ 2 
1296 9 



ln(a;) - -^{x + x'^)C2 ln(x) 



(8- 12x + 3x2 - 17x3)^2 



(A.6) 
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For the quark-gluon channel, we have: 



A^,A(^) = ^(2-2x + x^)ln^(l-x) 

(2592 - 2278,r- 111x2 - 288x'^) 2n (642 + 190 .r + 553 .r^) ^ 

■In [1 — X) m (1 — X) m{x) 



+ 



36 

(23887 - 17388 x - 2538 - 784 x^) 50 



18 



162 



^(2-2x + x2)C2 



ln(l - x) 



(1665 - 2040 x + 174x2 -400x3), ^ ^, ^ ^ 4(38 + 21 x + 39 x^) , ^ . 
+ ^ — ln(l - x) ln(x)+ ^ ln(l - x) ln2(x) 

- ^(46 + 298x + 139x2) ln(l - x) Li2(l - x) - 2(2 + 2x + x^) ln(l - x) Li2(l - x^) 

2..„ 2xx. X (302 + 474X + 339x2)^. / (1 - xV 
--(42- 142x-x2)Li3(l-x)-^ ^Lia' ^ 



X 



(2 + 2 



-4(2 + 2x + x' 



-Li3(l - x2 

, 1 — X 



(2 + 2,r + x2) 



Lis 



(l-x2) 



X^ 



1 + X 

(979 - 144x- 215x2 + 52x3 



1 — X 
1 + x 



+ 
+ 



18 

(166 + 222x + 33x2 + 4x'^ 
54 



T-n . ^ (142 + 374x + 245x2) 
-Li2(l - x) -I Li2(l - x) ln(x) 



Li2(l - x2) + 2(2 + 2x + x2) ln(x) Li2(l - x2 

2 _ 00«^3> 



(133 + 202 x + 115x2) ^ 3 (837- 1296 x + 234x2 -226x3) ^ 



ln3(x) 

27 ^ ^ 54 

(22042 - 32040 X - 5847x2 - 2464x3) (194 + 222x + 213x2; 



In^(x) 



<2 



ln(x) 



324 9 
(173719 - 156324 X - 12687x2 - 6148x3) (1071 - 710 x - 130x2 - 144x3) 



1944 



+ 



18 



C2 



+ ^(2-2x + x2)C3 



(A.7) 



and 



^U-)= ^^ ';/"'^ n2(l-x) 

- + 1,(1 _ ,) _ + 1,(1 _ (A.8) 

(2-2x + x2), 2. X (29 - 38x + 19x2), ^ ^ (265 - 418 x + 179x2) 
+ ^ ^ In^x) + ^ ln(x) + ^ 1 . 
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For the scattering of two identical quarks, we find 

KqAix) = -y (3 - 2x - a;2) ln2(l - x) - y (4 + 4a; + x^) ln2(l - x) ln(x) 

+ ^ (17 - 12 X - 5 x^) ln(l - x) + ^(12 - 8 X - 5 x^) ln(l - x) ln(x) 
o 9 

8 32 
+ -(4 + 4x + x^) ln(l - x) In^(x) - —{A + 4x + x^) ln(l - x)Li2(l - x) 

— (2-2x + x^)Li3(l-x)- — (50 + 46x + 13x2)Li3 ' ^ 

8 16 

- -(6 - 4x - x^)Li2(l - x) + —(19 + 17x + 5x2)Li2(l - x) ln(x) 

+ 1^(19 + 17 X + 5x2) - A (18 - 10 X - 9x2)ln2(x) 

2 16 

- ;r=(129 - 212 X - 69x2) ^^(^^^ _^ (4 + 43. + ^2^^^ ^^^^^^ 

+ 4(-86 + 53x + 33x2) ^ _ 2^ _ 



and 



9 ' ' ^ ' ^ 



+ ^(17 - 12x - 5x2)ln(l - x) + -(12 - 8x - 5x2)ln(l - x) ln(x) 
3 9 

8 32 

+ -(4 + 4x + x2) ln(l - x) ln2(x) - — (4 + 4x + x^) ln(l - x)Li2(l - x) 

9 9 



|(4 + 4x + x2)Li3 (-^) 



- -(6 - 4x - x2)Li2(l -x) + -(4 + 4x + x2)Li2(l - x) ln(x) 
9 o 

+ ^(4 + 4x + x2) ln3(x) - ^(3 - 2x - x^) ln2(x) 

2 16 

- -(43-68x-29x2)ln(x) + — (4 + 4x + x2)C2 ln(x) 
9 9 

2 32 

- -(61 - 46x - 15x2) + —(3 - 2x - x2) C2 



and 



AL^(x) = 0. (A.12) 



(A.9) 



For the scattering of a quark-(anti-)quark pair of distinct flavor, we find 

A^^,^(x) = - ^(3 - 2x - x2) ln2(l - x) - ^(4 + 4x + x2) ln2(l - x) ln(x) 

4 



(A.ll) 



55 
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For quark-anti-quark scattering (same quark flavor), we have 



„2 1 ooA ™3\ 



A^^^(x) = -^(14 + 21x-A8x^ + 13x^) ln^{l - x) - — (4 + 4a; + x^) ln2(l - x) ln(x) 
4 



+ — (75 + 892 X - 1351 x^ + 384 x"") ln(l - x) 
81 

+ — (76 + 72 X - 189 x^ + 64 x=^) ln(l - x) ln(x) 
81 

8 32 

+ -(4 + 4x + x^) ln(l - x) ln2(x) - — -(4 + 4x + x^) ln(l - x)Li2(l - x) 

9 9 

- ^(2 + 2x + x2)Li3(l - x) - ^(18 + 18x + 5x2)Li3 f-ii^ 
27 9 \ X 

+ ^(2 + 2x + x2)Li3(l - x2) + ^(2 + 2x + x2)Li3 (-^^^^ 



Li3 ( ^—^ ) - Lis ( -T— ^ 
1+x/ V l+x 



+ -(2 + 2x + x2) 

8 16 

(12 + 30 X - 93x2 - 26x^)Li2(l - x) + — (7 + 7x + 2x^)^2(1 - x) ln(x) 

81 9 

- (2x + x2 + 6x2)Li2(l - x2) - ^(2 + 2x + x2)Li2(l - x^) ln(x) 

+ 4; (5 + 5 X + x2) In^ (x) - J- (27 + 27 X - 81 x2 + 44 x^) ln2 (x) 
27 81 

2 16 
+ — (9- 1064 x + 2111x2 - 768x3)ln(x) + y(4 + 4x + x2)(2 ln(x) 

+ ^(101 - 462x + 520x2 - 159x2) + -^(11 + 93x - 147x2 + 43x3) C2 
81 81 



(A.13) 



and 



,p , , 32(1 -x)3 16(3 -9x + 12x2 -4x3) , ^ ^ 
= g-^ ^ ln(l - x) - ^ ln(x) 

_ 8(41 - 111 x + 111x2 -41x3) 
243 ■ 



(A.14) 



In the expressions above, we have identified the renorniaUzation and the factorization 
scale with the Higgs boson mass, /i/j = fip = Ma- The dependence on these scales is 
logarithmic and can readily be reconstructed by employing scale invariance of the total 
partonic cross section. 
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